UNCLASSIFIED 


AD  NUMBER 

AD005020 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  MAR  1953. 
Other  requests  shall  be  referred  to 
National  Aeronautics  and  Space 
Administration,  Washington,  DC. 

AUTHORITY 

NASA  TR  Server  Website 


THIS  PAGE  IS  UNCLASSIFIED 


- A**--*!* 


Reproduced  by 


firmed  Services  Technical  Information  flgencr 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING,  DAYTON,  2,  OHIO 


NACA  TN  2912 


NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS 

TECHNICAL  NOTE  2912 

THE  NORMAL  COMPONENT  OF  THE  INDUCED  VELOCITY  IN 
THE  VICINITY  OF  A LIFTING  ROTOR  AND  SOME 
EXAMPLES  OF  ITS  APPLICATION 
By  Walter  Castles,  Jr.  and  Jacob  Henri  De  Leeuw 
Georgia  Institute  of  Technology 


Washington 
March  1953 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
TECHNICAL  NOTE  2912 

THE  NORMAL  COMPONENT  OF  THE  INDUCED  VELOCITY  IN 
THE  VICINITY  OF  A LIFTING  ROTOR  AND  SOME 
EXAMPLES  OF  ITS  APPLICATION 
By  Walter  Castles , Jr.  and  Jacob  Henri  De  Leeuw 

SUMMARY 


This  paper  presents  a practical  method  for  computing  the  approximate 
values  of  the  normal  component  of  the  induced  velocity  at  points  in  the 
flow  field  of  a lifting  rotor.  Tables  and  graphs  of  the  relative  magni- 
tudes of  the  normal  component  of  the  induced  velocity  are  given  for 
selected  points  in  the  longitudinal  plane  of  symmetry  of  the  rotor  and 
on  the  lateral  rotor  axis . 

A method  is  also  presented  for  utilizing  the  tables  and  graphs  to 
determine  the  interference  induced  velocities  arising  from  the  second 
rotor  of  a tandem-  or  side-by-side-rotor  helicopter  and  the  induced  flow 
angle  at  a horizontal  tail  plane. 


INTRODUCTION 


This  work,  conducted  at  the  Georgia  Institute  of  Technology  State 
Engineering  Experiment  Station  under  the  sponsorship  and  with  the  finan- 
cial assistance  of  the  National  Advisory  Committee  for  Aeronautics,  was 
undertaken  in  an  attempt  to  obtain  a better  understanding  of  the  induced 
flow  in  the  vicinity  of  a lifting  rotor. 

Previous  investigations,  such  as  those  of  references  1 and  2,  demon- 
strated that  the  solution  of  the  integral  for  the  normal  component  of  the 
induced  velocity  at  the  center  of  the  rotor  could  be  obtained  in  an  ele- 
mentary form  provided  certain  approximations  were  made  as  to  the  distri- 
bution of  vorticity  in  the  wake.  However,  the  value  of  the  integral  for 
the  induced-velocity  component  at  an  arbitrary  point  in  the  rotor  flow 
field  cannot,  in  general,  be  expressed  in  terms  of  elementary  functions. 
Its  numerical  evaluation  for  a specific  case  presents  considerable 
difficulty. 

De  Leeuw,  in  reference  3;  investigated  the  feasibility  of  calcu- 
lating the  induced  velocity  at  arbitrary  points  in  the  vicinity  of  the 
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rotor  by  an  alternative  method  which  consisted  of  (l)  numerically  inte- 
grating the  increments  induced  by  the  vortex  ring  wake  elements  within 
a given  distance  of  the  point  and  (2)  summing  up  the  effect  of  the 
remainder  of  the  wake  by  an  approximate  integral.  This  approach  is  quite 
general  in  that  it  can  be  applied  to  any  wake  which  can  be  approximated 
by  an  assembly  of  vortex  rings . It  was  found  that  the  method  afforded 
satisfactory  accuracy  with  the  expenditure  of  a reasonable  amount  of 
effort,  since  the  values  of  the  normal  induced-velocity  component  for 
the  isolated  rings  may  be  precomputed  and  tabulated  for  repeated  use. 

The  scope  of  the  present  paper  is  limited  principally  to  a consid- 
eration of  the  values  of  the  normal  component  of  the  induced  velocity 
at  points  in  the  longitudinal  plane  of  symmetry  and  within  the  region 
likely  to  be  occupied  by  the  second  rotor  of  a tandem-rotor  helicopter. 

In  addition,  the  values  of  the  normal  component  of  the  induced  velocity 
were  calculated  for  points  on  the  lateral  axis  of  the  tip -path  plane 
over  the  distance  of  interest  for  the  case  of  a helicopter  with  laterally 
disposed  rotors. 

In  view  of  the  present  lack  of  experimental  evidence  as  to  the 
actual  wake  distribution  of  vorticity,  the  calculations  for  the  present 
paper  were  based  on  the  same  assumptions  for  the  wake  shape  as  those 
found  in  references  1 and  2.  These  assumptions  were  that  the  wake  vortex 
distribution  consists  of  a straight  elliptic  cylinder  formed  by  a uni- 
form, continuous  distribution  of  vortex  rings  of  infinitesimal  strength, 
lying  in  planes  parallel  to  the  tip-path  plane  and  extending  downstream 
to  infinity. 


SYMBOLS 


&0  constant  term  in  Fourier  series  for  blade  flapping  angle  p 

coefficient  of  cosine  term  of  Fourier  series  for  blade 
flapping  angle  p where 

P = aQ  - a^.  cos  - b1  sin  ^ - . . . 
b^  coefficient  of  sine  component  of  flapping  angle 

Crp  thrust  coefficient,  f/p 


drag  of  fuselage 
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nondimensional  shortest  distance  from  a point  P to  a 
vortex  ring;  \p  + (x  - l)2  (fig.  1) 

dg  nondimensional  largest  distance  from  a point  P to  a vortex 

ring,  \Jz2  + (x  + l)2  (fig.  1) 

E(t)  complete  elliptic  integral  of  first  kind 

K(t)  complete  elliptic  integral  of  second  kind 

R radius  of  vortex  ring;  also  radius  of  rotor 

Rp  radial  distance  of  a point  P from  axis  of  a vortex  ring 

(fig.  1) 

r nondimensional  radius  vector  in  rotor  XY -plane 

T rotor  thrust 

V velocity  of  helicopter  along  flight  path 

Vi  normal  component  of  velocity  induced  at  a point  P by 

whole  wake 

AVi  increment  of  normal  component  of  velocity  induced  at  a 

point  P by  that  portion  of  wake  which  is  beyond  the 
range  of  table  I 

v normal  component  of  induced  velocity  at  center  of  rotor 

vr  radial  component  of  velocity  induced  at  a point  P by  a 

vortex  ring 

vz  axial  component  of  velocity  induced  at  a point  P by  a 

vortex  ring 

W gross  weight  of  helicopter 

w slope  of  longitudinal  variation  of  nondimensional  induced 

velocity  in  plane  of  rotor 

X,Y,Z  rotor  axes  (fig.  2) 

nondimensional  radial  distance  of  a point  P from  axis  of 
a vortex  ring,  Rp^R  (fig*  1) 


x 
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x',y',z' 

nondimens ional  coordinates  of  a point  P with  respect  to 
rotor  axes  (fig.  3) 

y 

slope  of  lateral  variation  of  nondimens ional  induced 
velocity  in  plane  of  rotor 

Zp 

distance  of  a point  P from  plane  of  a vortex  ring 

Z 

1 

nondimens ional  distance  of  a point  P from  plane  of  a 
vortex  ring,  positive  in  direction  of  vz,  Zp/R 

(fig.  1) 

a 

angle  of  attack  of  plane  of  zero  feathering 

ai 

induced  angle  of  attack 

Op 

fuselage  angle  of  attack 

°v 

angle  of  attack  of  tip -path  plane 

Pp 

angle  between  radius  vector  from  center  of  rotor  to  a 
point  P lying  in  XZ-plane  and  positive  X-axis,  positive 
above  rotor  (fig.  2) 

r 

vortex  strength 

X = (V  sin  a - 

v)/ftR 

\v  = (V  sin  oty. 

- v)/nR 

M-  = V cos  a/ftR 

by  = V cos  Oy/aR 

P 

density  of  air 

T _ d2  “ dl 
d2  + dl 

K 

angle  between  flight  path  and  horizontal,  positive  below 
horizontal 

X 

angle  between  axis  of  the  wake  and  normal  to  tip-path 
plane  (fig.  2) 

""  1 
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ft 

Subscripts : 

B 

F 


azimuth  angle  measured  in  XY -plane  between  radius  vector 
to  a point  and  positive  X-axis,  positive  in  going  from 
positive  X-axis  to  positive  Y-axis 

angular  velocity  of  rotor,  radians/ sec 


values  of  back  rotor  of  two  rotors  in  tandem 

values  of  front  rotor  of  two  rotors  in  tandem 

values  taken  with  respect  to  virtual  axis  of  rotation  or 
to  tip-path  plane 


ANALYSIS 

Velocity  Induced  by  a Vortex  Ring 


It  is  shown  in  reference  4 (ch.  VII,  sec.  l6l,  p.  237)  that  the 
stream  function  at  a point  P (fig.  l)  in  the  flow  field  of  a vortex  ring 
of  strength  r and  radius  R may  be  expressed  as 

* = - —(dp  + cL2)[k(t)  - E(tJ|  (1) 

where  R is  the  radius  of  the  vortex  ring,  d^_R  and  dgR  are  the  least 
and  greatest  distances  of  the  point  P to  the  vortex  ring, 


T 


d2  ~ dl 
d2  + dl 


(2) 


and  K(t)  and  E(t)  are  the  complete  elliptic  integrals  of  the  first 
and  second  kinds,  respectively. 

The  flow  field  of  a vortex  ring  is  axially  symmetric  and  thus  the 
axial  and  radial  velocity  components  vz  and  vr  at  a point  P,  having 

an  axial  distance  Zp  from  the  plane  of  the  vortex  ring  and  a radial 

distance  Rp  from  the  axis  of  symmetry,  are  given  by 


1 m 

Rp  dRp 


(3) 
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and 


v = 

Rp  dZp 


w 


It  is  shown  in  reference  3 that  equations  (3)  and  (4)  may  be 
expressed  as 


v_  = -I—(AB  + CDF) 
z 2nxR  ' 


(5) 


and 


v = -^£-(AB ' + CDF  ’ ) 
r 2*xR 


(6) 


where 


> 

11 

H 

fd 

-1 

(T) 

B _ X - 1 rx  + l 

(8) 

dl  d2 

C = h + d2 

(9) 

d = MiL 

(10) 

1 - T 


_ (l  + x2  + z2)  - d]_d2  (l  + x)d^2  + (l  - x)d22 

2xdBd2 


2x£ 


(ID 


B'  = z 


(12) 


F*  = 


zL  1 + x2  + z2^ 
xl  " dxd2 


(13) 
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= \Jz2  + (x  - l)2 

(H) 

= Jz2  + (x  + 1)2 

(15) 

and 


x nondimens ional  radial  distance  of  P from  axis  of  vortex  ring,  Rp/R 

z nondimens ional  distance  of  P from  plane  of  vortex  ring,  taken  posi- 

tive in  direction  of  vz  on  ring  axis,  Zp/R 

The  values  of  vz  and  vr  given  by  equations  (5)  and  (6)  become 

indeterminate  for  points  on  the  vortex-ring  axis  where  x = 0.  In  this 
case  it  follows  from  the  symmetry  of  the  flow  that  the  radial  component 
of  induced  velocity  is  zero,  and  the  axial  component  of  induced  velocity 
is  shown  in  reference  5 to  be 


(16) 


Numerical  values  of  vzR^r,  which  is  a nondimens  ional  factor 
expressing  the  normal  component  of  the  Induced  velocity  vz  in  the 

vicinity  of  a vortex  ring,  are  given  in  table  I.  The  table  includes  a 
range  of  nondimens ional  axial  distances  of  -4.2  < z < 4.2  and  of  nondi- 

mensional  radial  distances  of  0 < x < 5.0.  The  increments  of  z at 
which  the  values  of  vzR^F  are  given  are  suitable  for  numerical  inte- 
gration by  Simpson's  rule.  The  tabulated  values  were  obtained  by  calcu- 
lation or  by  interpolation  as  indicated  in  the  table.  With  the  exception 
of  those  points  which  are  close  to  the  circumference  of  the  vortex  ring, 
the  calculated  values  are  accurate  to  four  places. 


Normal  Component  of  Induced  Velocity 
in  Vicinity  of  a Lifting  Rotor 

It  is  assumed  in  this  report,  as  in  references  1 and  2,  that  the 
rotor  wake  vortex  distribution  consists  of  a straight  elliptic  cylinder 
formed  by  a uniform  distribution  of  an  infinite  number  of  vortex  rings 
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of  infinitesimal  strength,  lying  in  planes  parallel  to  the  tip-path  plane 
and  extending  downstream  to  infinity.  The  above-described  vortex  distri- 
bution is  equivalent  to  a vortex  sheet  of  uniform  finite  strength  per 
unit  length  dr/dZ  measured  in  the  Z-direction.  This  sheet  forms  a 
straight  elliptic  cylinder  coinciding  with  the  boundary  of  the  wake . 

Within  the  limitations  of  the  initial  assumptions,  it  may  be  shown 
from  the  results  of  references  1 and  2 that 


dT 

dZ 


fiRCr 


(*  ■ 1 -!) 


(17) 


where  the  subscript  v denotes  values  with  respect  to  tip-path-plane 
coordinates . 

The  increment  of  the  normal  component  of  velocity  at  a point  P in 
the  vicinity  of  the  rotor,  induced  by  the  wake  vortex  rings  within  the 
distance  from  P covered  by  table  I,  may  thus  be  found  by  graphical  or 
numerical  integration.  This  increment  constitutes  about  95  percent  of 
the  total  value  of  the  normal  component  at  the  center  of  the  rotor  and 
a large  part  of  the  total  value  for  most  points  within  the  region  con- 
sidered in  this  paper. 

The  contribution  of  the  vortex  rings  beyond  the  range  of  table  I 
to  the  induced  velocity  at  P may  thus  be  summed,  with  small  error  in  the 
final  result,  by  an  approximate  expression  which  is  integrable. 

The  value  of  the  velocity  potential  Z$p  at  P due  to  a closed 

vortex  element  of  strength  T is  shown  in  reference  4 (ch.  VII,  sec.  150, 
p.  212)  to  be 


(18) 


where  o>  is  the  solid  angle  subtended  at  P by  the  closed  vortex  element. 

It  is  a good  approximation  for  those  wake  vortex  rings  at  distances 
from  P beyond  the  range  of  table  I that  the  subtended  solid  angle  at  P 
is  equal  to  three  times  that  volume  cut  off  the  cone,  determined  by  P 
and  the  ring,  by  a plane  which  is  parallel  to  the  plane  of  the  ring  and 
which  is  located  a unit  distance  from  P.  It  follows  that 
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’ 2T 


Mp  ~ i.  zxjR  dz 


Consequently,  the  increment  to  the  normal  component  of  velocity  induced 
at  P by  that  portion  of  the  wake  extending  from  the  limit  of  table  I, 
at  z = Zp  to  z = 00  may  be  obtained  from  the  integral 


AV<  = 1 P 
* R dz 


d 1 dr  zx2  dz 

52  5dZ(z2  + X2)3/2 


It  is  shown  in  reference  3 that  equation  (20)  may  be  integrated  to 
obtain  the  value  of  AV^  at  a point  P having  coordinates  x' , y', 

and  z'  from  the  center  of  the  rotor  and  that  the  result  is 


1 dT  /2'/e  2cz2  * b\/j 

2 4Z  U qjK  \ 


4ac  + b‘ 


(b2  - 2ac) z„  + ab 


where 


= (x1  - z'  tan  X)2  + (y')2 


b = -2x'  tan  X 


c = 1 + tan  X 


q = 4ac  - b* 


K = a + bz2  + cz^ 


Z2  = Z1  " z 


I 
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For  the  point  P(0,0,0)  the  value  of  AV^  given  by  equation  (21) 

becomes  indeterminate.  It  is  possible,  however,  to  substitute  the  zero 
coordinates  in  the  equation  before  integrating.  Doing  so  yields 

(AVi)o,0,0  = - jpjj  cos3X(3  cos2X  - l)  (28) 


The  normal  component  of  the  induced  velocity  at  any  point  P(x',y',z') 
may  thus  be  found  in  terms  of  dr/dZ  by  adding  the  increment  obtained 
from  the  numerical  integration  of  the  values  induced  by  the  wake  vortex 
rings  within  the  range  covered  by  table  I to  AV^  obtained  from  equa- 
tion (21)  or  (28). 

In  the  present  analysis  where  the  rotor  wake  vortex  distribution  is 
approximated  by  a straight  elliptic  cylinder  there  arises  the  question 
as  to  whether  the  wake  angle  should  be  taken  as  that  at  the  rotor  or  that 
in  the  ultimate  wake.  As  the  induced  velocity  distributions  in  the 
vicinity  of  the  rotor  are  more  sensitive  to  changes  in  position  of  the 
adjacent  vortex  elements  than  to  changes  in  position  of  the  vortex  ele- 
ments at  the  greater  distances,  the  wake  angle  at  the  rotor  will  be  used 
in  the  present  analysis. 

It  follows  from  figure  2 that,  for  X < 90° , or 
= A,  cos  a-^  + n sin  a^  < 0, 


and,  for  X > 90°,  or  = \ cos  a-L  + \i  sin  a-j_  > 0, 


In  the  above  equations  a-j_  is  the  coefficient  of  the  cosine  term  of  the 
Fourier  series  for  the  blade  flapping  angle 

P = Sq  - a^  cos  t - b^  sin  t - • • • 
where  p is  measured  from  the  plane  of  zero  feathering. 


fr 
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For  x = 90° ) equation  (21)  is  indeterminant.  However,  by  replacing 

dr/dz  in  equation  (20)  by  its  equivalent  (—)  tan  X after  performing 

\dx/ 

the  indicated  differentiation  with  respect  to  z,  it  can  be  shown  that 
for  this  wake  angle 


(31) 

where  the  integral  now  covers  the  region  from  x = x^  to  x = » and 

^p  ^RCrp 

^ My(l  - | iV*)  »(l  * f *2) 

RESULTS 

The  results  are  presented  in  the  form  of  tables  and  graphs  of  the 
ratio  of  the  normal  component  of  the  induced  velocity  Vp  at  any  point 

P(pp,x/R  or  (Jp,  Y/R),  as  in  figure  2,  to  the  normal  component  of  the 

induced  velocity  v at  the  center  of  the  rotor.  It  is  shown  in 
reference  2 that 


(32) 


Consequently,  the  value  of  at  P may  be  easily  computed  from  the 

values  of  Vp/v  in  the  tables  and  graphs. 

Table  II  gives  the  values  of  V^/v  for  -3*2  % X/R  ^3*2  and 
tan  Pp  = -l/2,  -l/4,  0,  l/4,  and  l/2.  Table  III  gives  the  values  of 

Vi/v  along  the  lateral  axis  of  the  tip-path  plane.  Figures  4(a)  to  4(i) 
show  the  lines  of  constant  values  of  V^/v  in  the  longitudinal  plane  of 
symmetry  for  wake  angles  having  tangents  of  0,  l/4,  l/2,  1,  2,  4,  00,  -4, 


i 
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and  -2.  Figures  5 and  6 show  the  variation  of  V^/v  with  X for  points 
on  the  longitudinal  and  lateral  axes  of  the  tip-path  plane. 


APPLICATION  OF  RESULTS 

Determination  of  Mean  Value  of  Normal  Component  of 
Induced  Velocity  Over  Front  and  Back  Rotors 
of  a Tandem-Rotor  Helicopter 


Making  the  approximation  that  the  mean  values  of  the  induced  velocity 
are  the  values  at  the  centers  of  the  respective  rotors,  and  being  given 
the  flight-path  velocity,  climb  angle,  gross  weight,  fuselage  drag, 
fuselage  angle  of  attack,  thrust  and  tip  speed  of  the  front  and  rear 
rotors,  and  the  geometry  of  the  helicopter,  the  mean  values  of  the 
induced  velocities  may  be  found  as  follows: 

The  angles  of  attack  ay 
rear  rotors  are  very  nearly 

Oy  — 

w - Uf  sm  fJc 


of  the  tip -path  planes  of  the  front  and 


Df*  cos  j6 


~r 


(33) 


where 

fic  angle  between  flight  path  and  horizontal,  positive  below  horizontal 
drag  of  fuselage 
W gross  weight 

Denote*  values  of  the  parameters  of  the  front  rotor  by  the  subscript  F 
and  of  the  back  rotor  by  the  subscript  B.  Then 


and 


% = (s)r  cos 


(3*0 


COS  Oy 


(35) 
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As  a first  approximation,  the  interference  induced  velocity  at  the  front 
rotor  due  to  the  thrust  of  the  back  rotor  may  be  neglected.  Then 


/V  \ VF 

( — sin  cl. 

\flR/F  (nR)p 


(36) 


where  the  value  of  vp 


is  given  by  equation  (32)  or,  for 


> 0.15, 


v 

HR 


(37) 


The  value  of  Xp  may  then  be  obtained  from  equation  (29)  or  (30) . When 
X,  o<y,  Of,  Vp,  and  the  geometry  of  the  helicopter  are  known,  the 


position  of  the  center  of  the  rear  rotor  with  respect  to  the  front  rotor 
can  be  determined.  Then  the  nondimens ional  velocity  V-j_/v  induced  at 


the 

may 

the 


center  of  the  rear  rotor,  because  of  the  thrust  of  the  front  rotor, 
be  obtained  by  interpolation  from  one  of  figures  4(a)  to  4(i)  for 


appropriate  value  of  Xp.  The  approximate  value  of  Vj^  ^ 


is  then 


vBtotal  ’ (VvW  + tb 


(38) 


where  vfi 
values  of 


can  be  obtained  from  equation  (32)  or  (37)*  The  approximate 

Xy  , X^,  and,  thus,  the  interference  induced  velocity  at  the 

B 


front  rotor  may  be  found  to  evaluate  v™  . In  general,  it  will  be 

rtotal 

necessary  to  iterate  v-d  for  an  accurate  result,  because  of  the 

rapid  variation  of  the  interference  induced  velocity  at  the  rear  rotor 
with  change  in  the  wake  angle  of  the  front  rotor  and  with  position  of 
the  rear  rotor  with  respect  to  the  tip-path  plane  of  the  front  rotor. 


For  a tandem-rotor  helicopter,  having  approximately  equally  loaded 
rotors  of  equal  size  spaced  approximately  1 rotor  diameter  apart  with 
small  vertical  offset  and  operating  in  the  high-speed  flight  range,  it 
is  seen  from  figures  4(f)  and  4(g)  that 
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and 


V*total  ^ °>^CT 

aE 


total 

or 


1.25Ct 


U, 


MO 


(39) 


m 


Determination  of  Longitudinal  Variation  of  Normal 
Component  of  Induced  Velocity  Over  Front  and 
Back  Rotors  of  a Tandem-Rotor  Helicopter 

If  the  normal  component  Vi  of  the  induced  velocity  at  P(r,V)  on 
the  rotor  disk  is  approximated  by  the  expression 


V, 

— = - — + yr  sin  \|r  + wr  cos  (hi) 

fiR  DR  J 


it  may  be  shown  from  the  results  given  in  reference  2 that,  for  a single 
rotor, 


The  increments  in  w arising  from  the  second  rotor  of  a tandem-rotor 
helicopter,  to  be  added  to  the  values  given  by  equation  (h2)  for  the 
front  and  rear  rotors,  may  be  obtained  in  the  general  case  from  the 
values  of  V-^/v  from  the  figures  at  r = 0.75  and  + = 0 and  « on 

the  respective  rotors.  Thus, 


GL  -G) 


^=0 
r=0.75 


V=jt 
r=0.75| 


Aw  = 


2 v_ 

3 fiR 


(43) 
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For  high-speed  flight  and  small  overlaps  between  the  rotor  disks 


Aw_,  « 7 — (44) 

F 6 HR 

and,  with  less  accuracy, 

AwB  ~ f — (45) 

B 4 fiR  ' 


Determination  of  Induced  Flow  Angle 
at  a Horizontal  Tail  Plane 


When  the  values  of  a^,  v/fiR,  X,  and  for  the  rotor  or  rotors 

in  question  have  been  determined  and  when  the  helicopter  flight  condi- 
tion under  consideration  is  known,  the  geometric  position,  and  thus  the 
values  of  Z/R  and  X/R,  of  the  horizontal  tail  plane  may  be  calculated 
and  the  value  or  values  of  V^v,  found  from  the  figures . Then  the 

induced  angle  cuj_  . at  the  tail  plane  is  approximately 
tail 


a 


Hail 


V cos  aj. 


(46) 


CONCLUDING  DISCUSSION 


The  assumption  that  the  planes  of  the  wake  vortex  rings  remain 
parallel  to  the  tip -path  plane  is  the  only  one  of  the  various  initial 
assumptions  as  to  the  wake  distribution  of  vorticity  which  appears  likely 
to  affect  the  engineering  accuracy  of  these  results  at  the  higher  flight 
speeds.  It  is  the  opinion  of  the  senior  author  that  the  present  investi- 
gation and  that  of  reference  1 indicate  that  the  planes  of  the  wake  vor- 
tex rings  must  be  tilted  to  the  rear  as  they  leave  the  rotor,  possibly 
approaching  a tilt  angle  in  the  ultimate  wake  of  half  the  wake  angle  X. 
The  quantitative  effects  of  such  a tilting  of  the  wake  vortex  rings  may 
not  be  large,  as  the  increments  of  the  radial  components  of  the  induced 
velocity  introduced  because  of  the  tilt  of  the  wake  vortex  rings  will 
tend  to  compensate  for  the  decrease  in  the  normal  component . 


\ 


» 


1 6 
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For  the  lower-speed  flight  conditions  the  initial  assumptions  as  to 
the  wake  distribution  of  vorticity  are  compatible  only  with  the'  assump- 
tion that  the  generating  rotor  is  lightly  loaded  and  has  blades  with 
constant  circulation  along  the  radii.  Caution  should  therefore  be  exer- 
cised in  applying  the  results  of  this  analysis  to  points  on  or  close  to 
the  disk  of  a specific  rotor  which  is  operating  at  the  lower  flight-path 
velocities. 

At  the  center  of  the  rotor,  where  the  values  of  the  induced  velocity 
calculated  by  the  method  presented  in  this  report  could  be  compared  with 
the  values  obtained  from  the  exact  integral,  the  error  was  in  every  case 
less  than  1 percent.  However,  for  those  points  in  the  flow  field  that 
lie  close  to  the  wake  vortex  sheet,  there  are  irregularities  in  the 
tabulated  values  of  V^/v  caused  by  difficulties  with  the  interpolations 

in  the  table  of  the  values  of  the  normal  component  of  the  induced  velocity 
for  the  vortex  rings . 

It  should  be  noted  that  the  present  analysis  neglects  the  effects 
of  the  lateral  dissymmetry  of  the  blade-bound  vortices  and  the  consequent 
lateral  dissymmetry  of  the  wake  vortex  elements  which  occur  in  forward 
flight.  Numerical  calculations  show  that  the  effects  of  these  lateral 
dissymmetries  on  the  induced  velocity  distributions  are  small  at  points 
of  interest  outside  the  boundaries  of  the  rotor  disk  or  wake  but  should 
be  taken  into  account  in  computing  the  longitudinal  and  lateral  distribu- 
tions across  the  rotor  disk.  The  first-order  effects  of  the  lateral 
variation  in  the  strength  of  the  blade -bound  vortices  are  accounted  for 
by  equation  ( 42 ) and  the  first-order  effects  of  the  lateral  variation 
in  the  strength  of  the  wake  vortices  may  be  taken  into  account  by  use  of 
a value  of 


y 


2h 


v 

fiR 


in  equation  (4l) . 

In  order  to  construct  figures  4 to  6 it  was  necessary  to  compute  a 
large  number  of  values  of  V^v  in  addition  to  those  listed  in  table  II. 

However,  as  these  additional  points  were  at  scattered  locations,  and 
consequently  of  little  use  for  any  other  purpose,  they  have  been  omitted 
from  this  report. 

Since  the  wake  angles  of  the  rotors  of  helicopters  operating  in  the 
upper  half  of  their  speed  range  fall  in  a narrow  band  between  80°  and  85°, 
it  would  be  useful  to  have  the  induced  velocity  distribution  for  a wake 
angle  of,  say,  82°.  However,  investigation  showed  that  in  order  to 
obtain  sufficiently  accurate  values  for  this  wake  angle  it  would  first 


MCA  TN  2912 


17 


be  necessary  to  compute  a large  number  of  additional  values  of  the 
induced  velocity  of  the  vortex  ring  for  the  region  within  two-tenths  of 
a ring  radius  from  the  periphery  of  the  ring.  The  computations  for  the 
82°  wake  angle  were  therefore  too  lengthy  for  the  results  to  be  included 
in  the  present  report. 

It  appears  from  the  results  of  this  investigation  that  the  inter- 
ference induced  velocity  at  the  rear  rotor  of  a tandem-rotor  helicopter 
in  high-speed  flight,  due  to  the  thrust  of  the  front  rotor,  is  of  the 
same  order  of  magnitude  and  of  the  same  sign  as  the  self -induced  velocity. 
Consequently,  the  interference  induced  velocity  should  be  taken  into 
account  in  longitudinal  stability  calculations  and  in  computing  the 
equilibrium  values  of  the  mean  blade  angle  and  torque  coefficients.  The 
interference  induced  velocity  at  the  front  rotor  of  a tandem-rotor  heli- 
copter in  high-speed  flight  is  of  the  order  of  7 percent  of  the  self- 
induced  velocity  and  is  opposite  in  sign. 

The  longitudinal  gradient  of  the  interference  induced  velocities 
at  both  rotors  of  a tandem-rotor  helicopter  in  high-speed  flight  is  of 
.opposite  sign  to  the  longitudinal  gradients  of  the  self-induced  veloc- 
ities, and,  consequently,  will  have  the  effect  of  reducing  the  required 
equilibrium  values  of  b^,  the  coefficient  of  the  sine  component  of  the 

flapping  angle. 

For  side -by- side -rot or  helicopters  in  high-speed  flight  the  mean 
values  of  the  interference  induced  velocities  are  of  the  order  of  15  per- 
cent of  the  self-induced  velocities  and  are  opposite  in  sign.  The  lat- 
eral gradients  of  the  mutual  interference  induced  velocities  are  large 
for  the  adjacent  portions  of  the  rotors.  These  large  gradients  may 
cause  early  tip  stall  if  the  rotor  rotation  is  such  that  the  retreating 
blades  are  in  this  adjacent  rotor  position. 

The  normal  component  of  the  induced  velocity  inside  the  wake  of  a 
helicopter  rotor  in  high-speed  flight  appears  to  reach  its  maximum  and 
final  value  of  about  twice  the  value  at  the  center  of  the  rotor  at  a 
distance  of  about  1 rotor  radius  downstream  from  the  center  of  the  rotor. 
For  hovering  and  the  lower  forward  speeds  the  induced  velocity  inside 
the  wake  reaches  its  final  value  of  about  twice  that  at  the  center  of 
the  rotor  at  a distance  of  about  2 rotor  radii  downstream. 


Georgia  Institute  of  Technology 
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TABLE  II 


NONDIMENSIONAL  VALUES  OF  NORMAL  COMPONENT  OF  INDUCED  VELOCITY  IN  LONGITUDINAL 
PLANE  OF  SYMMETRY  OF  A LIFTING  ROTOR  FOR  X £ 9O0 
[For  flight  conditions  for  which  X > 90°  and  a developed  wake 


exists,  values  of  V^v  for  ^X,  X/R,  are  the  same 
as  those  for  ^180°  - X,  X/R,  -|3pj 


V,/v  for  values  of  tan  fo  of  - 

x/r 

' 

-1/2 

-1/4 

0 

1/4 

1/2 

X 

= 0° 

(a) 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.40 

1.220 

1.112 

1.000 

.888 

.780 

.80 

1.555 

1.368 

1.000 

.632 

.445 

.90 

1.00 

1-643 

1.495 

1.000 

.500 

.505 

.363 

.357 

.282 

1.10 

-.224 

-.253 

0 

.253 

.224 

1.20 

-.180 

-.180 

0 

.180 

.180 

1.60 

-.088 

-.067 

0 

.067 

-.088 

2.00 

-.053 

-.038 

0 

.038 

.053 

3-20 

-.019 

-.012 

0 

.012 

.019 

Vj_/v  for  values  of  tan 

1 

<H 

O 

x/r 

-1/2 

-1/4 

0 

1/4 

1/2 

X = 14.04°  (tan  X = l/4) 

(b) 

-3.20 

0.011 

0.001 

-0.010 

-0.020 

-0.023 

-2.00 

.028 

.012 

-.027 

-,06l 

-.063 

-1.60 

.057 

.026 

-.046 

-.101 

-.105 

-1.20 

.132 

.099 

-.108 

-.232 

-.200 

-.80 

.363 

.521 

.860 

1.245 

1.406 

-.40 

• 729 

.835 

.948 

1.059 

1.170 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.40 

1.262 

l.l6l 

1.052 

.940 

.828 

.80 

1.626 

1.479 

1.139 

.746 

.526 

1.20 

1.40 

-.042 

-.036 

.150 

.281 

.242 

1.60 

.004 

■079 

.121 

* .126 

2.00 

-.013 

.007 

.046 

.073 

•077 

3.20 

-.001 

.004 

.016 

.034 

.028 

aValue  of  V-^/v  changes  2.000  units  in  passing  through  boundary  of  wake.  Values 

are  axially  symmetric  about  rotor  axis  and  antisymmetric  about  tip-path  plane;  values 
about  which  they  are  antisymmetric  are  1.000  for  x/R  < 1 and  0 for  X/R  > 1. 
bValue  of  V^/v  changes  I.9I+O  units  in  passing  through  boundary  of  wake. 


x/r 


-1/2 


lA 


1/2 


Vj/v  for  values  of  tan  fp  of  - 


-lA  0 

X = 26.56°  (tan  X ---  l/2) 
(c) 


-3-20 

0.005 

-0.004 

-0.016 

-0.024 

-0.026 

-2.00 

.019 

-.006 

-.044 

-.067 

-.065 

-1.60 

.034 

-.016 

-.073 

-.120 

-.109 

-1.20 

.089 

.038 

-.180 

-.261 

-.198 

-.80 

• 305 

.426 

•738 

-.4o 

.687 

.788 

.899 

1.046 

1.126 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.40 

1.314 

1.211 

1.101 

.987 

.874 

.80 

I.692 

1.582 

1.261 

.860 

.604 

1.20 

1.803 

.328 

.403 

• 305 

1.60 

.113 

.125 

.167 

.191 

.166 

2.00 

.088 

.084 

.105 

.115 

.105 

3-20 

.038 

•037 

.041 

.044 

.040 

V./v  for  values  of  tan  pp  of  - 

v / 

-1/2 

-1/4 

0 

1/4 

1/2 

X = 45.00°  (tan  X = l) 

(a) 

-3.20 

-0.002 

-0.010 

-0.019 

-0.025 

-0.023 

-2.00 

-.001 

-.028 

-.060 

-.071 

-.060 

-1.60 

.006 

-•039 

-.107 

-.122 

-.096 

-1.20 

.036 

-•039 

-.265 

-.249 

-.172 

-.80 

.217 

.285 

.535 

-.142 

-.40 

.622 

.719 

.824 

■ 927 

1.046 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.4o 

1-381 

1.283 

1.176 

.1.062 

.947 

.80 

1-782 

1.711 

1.465 

1.048 

.740 

1.20 

1-957 

1.866 

.683 

.593 

.423 

1.60 

1.829 

.467 

.410 

• 332 

.247 

2.00 

• 372 

.272 

.215 

.164 

3-20 

•307 

.165 

• 113 

.087 

.066 

Vj/v  for  values  of  tan  pp  of  - 

x/R 

-1/2 

-1/4 

0 

1/4 

1/2 

X = 63.430  (tan  X = 2) 
(e) 


-3-20 

-0.008 

-0.021 

-0.022 

-0.020 

-0.020 

-2.00 

-.015 

-.043 

-.064 

-.069 

-.051 

-1.60 

-.020 

-.067 

-.119 

-.116 

-.083 

-1.20 

-.012 

-.105 

-.294 

-.238 

-.140 

-.80 

.124 

.140 

.326 

-.320 

-.144 

-.40 

.550 

.639 

•746 

.858 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.40 

i.46o 

1.366 

1.254 

1.144 

1.032 

.80 

1.886 

1.849 

1.674 

1.265 

.911 

1.20 

2.022 

2.098 

1.240 

.900 

.600 

1.60 

2.030 

2.098 

.878 

.582 

■ 384 

2.00 

2.026 

.653 

.401 

.262 

3-20 

2.017 

•776 

.278 

.181 

.100 

cValue  of  V^v 
dValue  of  V^v 


changes  1.789  units  in  passing  through  boundary  of  vake. 
changes  1.4l4  unitB  in  passing  through  boundary  of  vake. 


'Value  of 


Vj^v  changes  0.894  units  in  passing  through  boundary  of  vake. 


Vp/v  for  values  of  tan  Pp  of  - 

X/  R 

-1/2 

-1/4 

0 

1/k 

1/2 

x = 75. 97' 

D (tan  X = 4) 

(f) 

-3.20 

-0.011 

-0.019 

• -0.026 

-0.022 

-0.014 

-2.00 

-.025 

-.052 

-.072 

-.065 

-.043 

-1.60 

-.035 

-.084 

-.124 

-.101 

-.067 

-1.20 

-.043 

-.136 

-.309 

-.213 

-.110 

-.80 

.062 

.054 

.166 

-.193 

-.087 

-AO 

.485 

.549  ! 

.675 

.48  3 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

AO 

1.494 

1.427 

1.325 

1.165 

1.096 

.80 

1.576 

1.942 

1.834 

1.436 

1.053 

1.20 

1.895 

2.150 

1.183 

.763 

1.60 

1.250 

2.086 

1.384 

.869 

.530 

2.00 

.967 

2.063 

1.187 

.623 

.377 

3.20 

.735 

2.022 

.785 

• 328 

.169 

V±/v  for  values  of  tan  Pp  of  - 


-1/4 

0 

1/4 

1/2 

x = 90.00° 

(g) 


-3.20 
-2.00 
-1.60 
-1.20 
i no 

-0.016 

-.035 

-.052 

-.076 

-0.022 

-.060 

-.097 

-.182 

-0.026 

-.078 

-.137 

-.331 

- 6?l 

-0.022 

-.060 

-.057 

-.182 

-0.016 

-.035 

-.052 

-.076 

-.80 

-.010 

-.057 

-.075 

-.057 

-.010 

-.40 

Al4 

.519 

.564 

.519 

.414 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

.40 

1.193 

1.312 

IA36 

1.312 

1.193 

.80 

1.266 

1.664 

2.075 

1.664 

1.266 

1.20 

1.046 

1.602 

2.331 

1.602 

1.046 

1.60 

.804 

1.354 

2.137 

1.354 

.804 

2.00 

.619 

1.175 

2.078 

1.175 

.619 

3-20 

.318 

•774 

2.026 

• 774 

.318 
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VORTEX-RING  AXIS 
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Figure  1.-  Coordinates  for  vortex  ring  and  table  I 
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Figure  4.-  Lines  of  constant  values  of  isoinduced  velocity  ratio 

in  longitudinal  plane  of  symmetry. 
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Figure  4.-  Continued. 
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